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Abstract 

We consider r-lepton mass effects in the cascade decays H —)• Z{-^ + Z*(—>■ 

T+r“) and H —)• W~{-^ Since the scale of the problem is 

set by the off-shellness of the respective gauge bosons in the limits (m^ + m^/)^ < 
< {niH — rnw,z)‘^ and not by lepton-mass effects are non-negligible for the 

r modes in particular close to the threshold of the off-shell decays. Lepton-mass 
effects show up in the rate and in the three-fold joint angular decay distribution for 
the decays. Nonzero lepton masses lead to leptonic helicity-flip contributions which 
in turn can generate novel angular dependencies in the respective three-fold angular 
decay distributions. Lepton-mass effects are more pronounced in the H —)• Z{-^ 
ii)Z*{-^ tt) mode which, in part, is due to the fact that the ratio of lepton helicity- 
flip/nonflip contributions in the decay Z* is four times larger than in the 

decay W~^* —>■ Overall the inclusion of r mass effects leads to a 3.97% reduction 

in the leptonic H ^ ZZ* rate. Lepton mass effects are quite pronounced for (f‘ values 


from threshold up to ~ 200 GeV^. For example, at = 50GeV^ the transverse- 
longitudinal-scalar helicity composition of the off-shell Z-boson changes from 0.06 : 
0.94 : 0 to 0.04 : 0.65 : 0.31 for the r lepton. This has observational consequences 
for the angular decay distributions of the final-state leptons. We also briefly consider 
the corresponding off-shell - off-shell decays H —)■ Z*{-^ -|- Z*(—)■ r+r“) and 

H T+Vr). 
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1 Introduction 


We consider lepton-mass effects in the off-shell decays of gange bosons in the processes 
Z* —>■ r+r“ and W^* —?• r+z/^ where the off-shell gange bosons are prodnced 

in the Higgs decays H —)• ZZ*, In the H ^ ZZ* case the corresponding 

£ = e, fi modes have recently been observed at the LHC and are therefore adeqnately 
dnbbed “Higgs discovery channels” BEI. Fnrther evidence on these decays has been 
presented in Ref. |3]. The qnantnm nnmbers of the Higgs boson have been pinned down 
by an angnlar analysis of the fonr leptons in the hnal state to be = O’*" both in the 
leptonic H —)■ ZZ* mode [31 HJ [5] as well as in the leptonic H —)■ W~W~^* mode [6]. On 
the theoretical side there have been a nnmber of papers analyzing the qnantnm nnmbers 
of the Higgs boson throngh an angnlar analysis of the fonr-lepton hnal state among which 
are Refs. [ziEiEiiiniEiiiaiisiEiiisiiniiniiiHicn]. The physics of the Higgs boson in 
all its aspects has been nicely reviewed in three recent papers [20112IIE2]. 

Off-shell effects in the decays involving massive leptons will lead to additional scalar 
and scalar-longitndinal interference contribntions well familiar from nentron beta decay, 
the semileptonic decay [23], or from the decays B —)■ [231 ISSj 

and Af, —)• AcTUt- [26|. The scalar and scalar-longitndinal interference contribntions are 
qnadratic in the lepton masses and can thns be neglected at the scale rri^z- However, for 
the off-shell decays H ^ ZZ*, W~W~^* the scale is not set by ^ bnt by the off-shellness 
of the respective gange bosons which extends from threshold (maximal 

recoil point) to the zero recoil point at = {mn — mw,zY, i-e. one has 

{me + mef<q^<{mH-rnw,zY- ( 1 ) 

One will therefore have to carefnlly consider r-lepton mass effects particnlarly in the g^ 
region close to threshold given by g^ = 4m^ and g^ = for the leptonic modes in 
the decays H ^ ZZ* and H —>■ WW*, respectively. Lepton-mass effects rednce the 
overall rate relative to the zero lepton-mass case. In addition, lepton-mass effects lead to 
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leptonic helicity-flip contributions which in turn can generate novel angular dependencies 
in the respective angular decay distributions. These angular dependencies can mimic new 
angular terms introduced by higher dimension effective coupling terms [HIlslIIE] or non- 
SM (HVV) coupling terms [TTl HH HI], r-lepton mass effects should therefore not be 
neglected if one is aiming for high precision physics in the Higgs sectorll] 

Our paper is structured as follows. After this introductory section, in Sec. 2 we present 
a general formula for the three-fold angular angular decay distribution for the on-shell - off- 
shell decays H —)■ VV* —)■ ££££. The angular decay distribution is obtained using helicity 
methods. In Sec. 3 we discuss lepton-mass effects in the decay H ^ ZZ* ^ ££tt and their 
effect on the rates and the angular decay distributions. We do the same in Sec. 4 for the 
decays H —)■ W~W~^* —)■ £v(,TVr- In Sec. 5 we summarize our results and conclude with 
some general remarks. Some technical material regarding helicity amplitudes is relegated 
to the Appendices. In Appendix A we list the helicity amplitudes for the H —)• VV* 
transitions. The helicity representation of the lepton tensors in the neutral- and charged- 
current cases can be found in Appendix B and C, respectively. 


2 General formalism 


The three-fold angular decay distribution in the cascade decays H —)■ VV* —)■ ££££, V = 
Z, W can be derived from the covariant contraction of the on-shell and off-shell lepton 
tensors and with the (HVV) Higgs coupling Ha/s where the vertices are connected 
by the propagator projectors (spin 1) and Pq^i (spin 0 © 1). One has 


W( 0 „ e„ X) = prip) PSm(q) r?j(p) LfAi) Pi^ip) Pml(i) 


0 / 3 ' 


( 2 ) 


where, in the Standard Model (SM), Haa' = Qaa'- We denote the on-shell and off-shell mo¬ 
menta of the gauge bosons by p and q. In the unitary gauge the on-shell spin-1 propagator 
^Lepton-mass effects in the rate H ^W£v and H —>■ ZU are also taken into account in Ref. m- 
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= iTT'v) and the off-shell propagator {q^ 7 ^ rriy) read 


pr{p) = -9^^+ 


pO■p^l 


P&(q) = -s"' + Pr- 


( 3 ) 


rriy 

Note that in the nnitary gangqj the off-shell propagator Po^iiq) contains a spin-1 and a 
spin-0 piece. This can be seen by splitting the off-shell gange propagator in Eq. ([3]) into 
its spin -1 and spin -0 components according to 


PSiiM = -f" + 


ui3 . q q 


m 


V 


= -9^^ + 


iy/3 , ^ ^ 


qPq!^ 


Fs{q% 


( 4 ) 


spin 1 


spin 0 


where 

Psif) = (1 - 1 ^) ■ ( 5 ) 

In the zero lepton-mass approximation one has = 0 and therefore the spin-0 piece 

in Eq. (|1]) does not contribnte and can be dropped when evalnating Eq. (|2]). This is 
always a good approximation for i = e, fi bnt no longer a good approximation for i = t. 
An interesting observation concerns the spin-0 contribntion. Taken together with the 
propagator pole proportional to — my)~^, the contribntion of the spin -0 piece can be 
seen to be proportional to a contact interaction of the form {HV'ijj'ip) with a g^-dependent 
conpling when one sets T^ = 0 . 

Technically there are two rontes to obtain angnlar decay distribntions from Eq. ([2]). In 
the hrst ronte one parametrizes the fonr-vectors of the problem in terms of the hve phase- 
space variables = rriy, q^ , cos 6p, cos 6q and y (cf. Figs. [Hand[7]). The covariant evalnation 
of the Lorentz-invariant expression ([2]) leads to a nnmber of scalar prodncts of momenta 
that are dehned in different reference frames. When doing the reqnisite contractions, the 
^The choice of the unitary gauge is mandatory to obtain a gauge-independent result. This can be seen by 
considering a general covariant gauge where one has to consider Goldstone boson exchange in addition 
to gauge boson exchange. In the coupling to the final state fermion pair the gauge parameter ^ cancels 
between the Goldstone and gauge boson contributions, resulting in the unitary gauge propagator. This 
has been explicitly demonstrated for fermion-fermion scattering and for the decay t ^ h + W'^* [55]. 
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four-momenta have to be boosted to a common reference frame as e.g. described in Ref. [13] 
for the decay H —)■ Zii and in Ref. [301 E] for fhe decay —)■ 7r^7r‘’e’''e . One then 
arrives at the desired three-fold joint angular decay distribution. 

A second, perhaps more intelligent route, is to use an analysis in terms of helicity am¬ 
plitudes. The advantage of the helicity method is that the origin of the angular factors 
multiplying the helicity structure functions can be straightforwardly identihed. The angu¬ 
lar factors can be seen to arise from the transformation properties of the helicity amplitudes 
under the action of the rotation group. 

In order to transform to the helicity representation of the covariant form in Eq. (|2!) 
one makes use of the completeness relation for the spin -1 on-shell and off-shell polarization 
vectors. The on-shell and off-shell propagator can be expanded according to [^ l32] 

pr(p) =-9“"+^ = i: E”(Av)f'‘(Ao (6) 


= my) and 


I f L-[ 


^ {Xv*)e*°‘{Xv*) gxy^Xy,- 


( 7 ) 


w Av*=i,±l,0 

Note that there is an additional spin-0 degree of freedom propagating in the off-shell 
propagator in Eq. ([7|). We shall specify this spin-0 degree of freedom by assigning the 
label Xv = t {t for time-component) to this mode. According to the separation in Eq. (|3|) 
the “t” mode carries the weight F 5 = (1 — /rriy) which hnally leads to g\y,\y, = 
diap {F 5 , —1, —1, —1} in Eq. ([7|). The four polarization four-vectors ±1,0) will be 

specihed in Appendix A. 

In low-energy calculations such as neutron /3 decay or in the semileptonic bottom- 
hadron decays one usually drops the term proportional to {q'^/my) in Eq. ([5|) since one has 
-C rriy. However, in the present application the factor can become as large as 

30% at the zero recoil point and can therefore not be neglectedjf 


^In muon decays fi —^ e + g-fi + Ve where one is aiming for ultrahigh precision, the importance of the 
contributions have been discussed in the literature 
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Figure 1: Definition of the momenta p and q, the polar angles 9p and 9q, and the azimuthal 
angle x in the cascade decay H —)■ e+e“) + Z*(—>■ 

With the help of the completeness relations (jd]) and ([7]) the covariant form of the angular 
decay distribution ([2]) can be cast into a representation in terms of helicity components. 
One has 


W{9p,9q,x) 


E i-^sY ^ ■^'L^ly^{cos9p 

Xy,\'^ = ±l,0 

( 8 ) 


)H 


Xy ,A^ 


.H 


\' 


' (g) 

X' 


(COS0g,X), 


where J = 1 for Ay = ±1,0, J = 0 for \y = t and correspondingly for the primed 
quantities. It turns out that J = J' in the decay H Z(—^ e+e“) + Z*(—>- r+r“) as long 
as one is not analyzing r-polarization effects, i.e. there are no spin -0 - spin -1 interference 
effects in this decay. 

The evaluation of the helicity components of the ih —)■ VV* transition amplitudes 
Hxy^Xy, is given in Appendix A while the evaluation of the helicity components of the 
lepton tensors (cos^p) and L^xl* x' (cos 6 'q,x) are given in Appendix B (neutral- 

current case) and C (charged-current case). 

Up to this point we have allowed for a general structure of the {HVV) coupling. In 
the following we shall specify to the SM coupling with Haa' = Qaa'- 
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3 The four-body decay H —> Z[-^ ) + Z*(—> r+r ) 


In this section we write down the three-fold angular decay distribution of the decay H —)■ 
Z(—)■ £+£“) -I- Z*(—)■ involving two different pairs of leptons, i.e. we assume i ^ t. 

The corresponding decay H Z{-^ £+£“) -|- Z*{-^ i^i~) involving two pairs of identical 
leptons (with and without lepton-mass effects) is more difficult to analyze due to the 
presence of nonfactorizing interference contributions. These identical-particle effects will 
be treated in a separate paper [35] . 

3.1 Three-fold angular decay distribution for the 
four-body decay H —> Z{^ + Z*(—)• r+r“) 

We begin our discussion by presenting an explicit form of the three-fold angular decay 
distribution given by Eq. ([8|). The relevant helicity components of the on-shell and off- 
shell lepton tensors are listed in Appendix B while the helicity components of the H ^ ZZ* 
transition amplitude can be found in Appendix A. The polar angles 9p and 6q are defined in 
the respective lepton pair center-of-mass systems as shown in Fig. [1] The azimuthal angle 
X describes the relative orientation of the two decay planes. We split the decay distribution 
into a helicity-nonflip and helicity-flip part, 

{2p^2q^)-^W^f{ep,eq,x) = {P++ + P-) 

X (|(1 + cos^0p)(n| -f a'jvp){l -f cos^ 6q){v'j + ajv'^) + AcosOpCosdgVga'jvpVq^ 

+ Poo sin^ Opivj + ajvl) sin^ 9q{vj + -F (p++ - p__) 

X (^(1 -I- cos^ 9p){vj -t- a^Up) cos9qVq + cos9pVp{l + cos^ 9q){vj + ajvq)'^ V£ae 
-h (p+o -h p-o) sin9p sin9q(4:V£ajvpVq + cos9p{vj -h a^v^) cos9q{vj + a^Vq)^ cosy 
-h 2(p+o - P-o) sin 9p sin 9qVeae (cos 9p{vj -1- a‘jvp)vq + cos 9q{v£ + ajvq)vp'^ cos y 
-h |p+_ sin^ 9p{vj + a]vp) sin^ 9q{vl + a^n|) cos2x 
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( 9 ) 


and 


(2p22g2) ^W^f{ep,eg,x) = -^{{v]+ a^ivl) X 

Q(p++ + p —)(1 + cos^ Op) sin^ Oqv'^ + poo sin^ Op cos^ Ogvj + ps Fg sin^ Opoj 
— \{p+o + p-o) sin 20p sin 20gv‘j cos x — \p+- sin^ Op sin^ Ogvj cos 2 %) 

+ (p++ — p _) cosOpVpSw? Ogvfai — (p+o — p-o) sm.0pVpSm.20gvlaicosx^, ( 10 ) 


where = 1 — 4mjp/p‘^ and = 1 — For symmetry reasons and for later 

applications in the off-shell - off-shell case we have written p^ for m| and = 1 — 4m|p/p^ 
for Wp = 1 on the on-shell side. The donble spin-density matrix elements Pmm' are bilinear 
forms of the helicity amplitndes describing the H ^ ZZ* transitions. They are dehned in 
Appendix A. 

In Eqs. dl]) and flTOj) we have also inclnded the contribntions from the parity-violating 

terms proportional to (p++ — p _) and (p+o — P-o)- These coefficient fnnctions are not 

popnlated by the parity-conserving SM [HVV) conpling. In Appendix A we briefly discnss 
the contribntion of a parity-violating non-SM conpling proportional to which 

wonld popnlate the (p++ — p_ ) and (p+o — P-o) coefficient fnnctions |T3, Ell [IH]. 

We add the flip and non-flip contribntions and expand the resnlt in terms of the Leg¬ 
endre polynomials Pi (cos 6 *) = cos 6 * and P 2 (cos 6 *) = |(3cos^6* — 1). The resnlt is written 
in the form 


(2p^2gX'W^{e„e„x) = ^E (/f + £9f)Ai(9p,9„,\), (11) 

^i=0 ^ “ 


i=0 


where in the second eqnation of ffTTl) we have split the coefficient fnnction Ff into its 
helicity-flip and helicity-nonflip part nsing the notation e = rn^/q^. 

The coefficient fnnctions and gf and their associated angnlar factors hi{0p,0g,x) 
are listed in Table [T] where we nse the abbreviation = yj -p ajvf with i = p,q. In 
addition we nse a short-hand notation for the donble density matrix elements, namely 
pu = P++ + p—, Pl = Poo, Pu+L = Pu + Pl and ps = ptt- Note that we have dropped the 
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i 

/f 

9i 

hi{9p^ 9q, y) 

0 

qSciS pu+L 

‘^C'£>{vjpu+L + ^ajPlps) 

1 

1 

^2C^Mljipu - 2 pl) 

-2Ci^vj{pu - 2pl) 

P2{C0S 9q) 

2 

^2C^Mljipu - 2 pl) 

-2pl)- QajFips) 

P2{cOS 9p) 

3 

iC<gCS(p„+4pi) 

-Ci^vj{pu + 4pi) 

P2{cOs9p)P2{cOS 9q) 

4 

QvjajvpVgPu 

0 

COS 9p cos 9q 

5 

9v‘^ajvpVq{p+o + p_o) 

0 

sin 9p sin 9g cos y 

6 

+ P-o) 

-|C'iS?^'l(P+o + p-o) 

sin 29p sin 29g cos y 

7 

|c<gc<»)p+_ 

-|c'iS?^'|p+- 

sin^ 9p sin^ 9g cos 2y 


Table 1: Coefficient functions appearing in the three-fold angular decay distribution of the 
decay H Z*%-^ £+£") + r+r") 

contributions of the parity-violating terms proportional to (p++ — p _) and (p+o — P-o) in 

Table m which are not populated by the parity-conserving SM (HVV) coupling. 

The dominant flip contributions proportional to aj are contained in and pf. Com¬ 
pared to ae the leptonic vector coupling Vi = -1-1- sin^ is much suppressed. This is 
different in the quark-antiquark case treated in Refs. [36l |37j where the electroweak vector 
and axial vector couplings to the quark pairs are comparable in size. As a result the pattern 
of the helicity-flip contributions in the quark pair production case is quite different from 
the lepton-pair production case [Ml [37] . 

In order to save space we have not expanded the angular factor sin^ 9p sin^ 9g in the last 
row of Table [1] The relevant expansion would be given by 

sin^ sin^ 6*, = ^ - Palcos^p) - P2{cos9g) + P2{cos9p)P2{cos9g)^. (12) 
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We then define a normalized decay distribution 


w^{e,,e,,x) = 




jw^{e'p,e'^,x')dcose'pdcose'^dx' sttV 

where Tf = J=f (and jf = /f/-F q^, = gf! J^q) and where 




X 


(13) 


+ eg^, = Ci^Ci^pu+L + 2eCi^iv^,pu+L + ^alFlps). 


(14) 


The normalized angular decay distribution {Op^Oq^x) obviously integrates to 1, i.e. 

J {9p,9q,x) d cos 9pd cos 9qdx = 1- (15) 

Before we start discussing our numerical results we want to specify our mass, width 
and coupling input parameters. We use the central value of the Higgs mass mn = 
125.09(24) GeV from the combined ATLAS and CMS measurement [3H]. For the remaining 
parameters we use the central values from the PDG [59] given by 


mw = 80.385(15) GeV, Tw = 2.085(42) GeV, 
mz = 91.1876(21) GeV, Tz = 2.4952(23) GeV, 

TUr = 1.77682(16) GeV, 

sin^ 9w = 0.23126(5), Gp = 1.1663787(6) x lO"® GeV^. (16) 

Our formulas are written in terms of the dimensionless coupling constant g^ which is related 
to Gp hj g‘^ = 8m^Gp/\/2. For practical numerical purposes we choose rriip = irie (or 
= 0) on the p side. On the off-shell q side we write meg = which can take the values 
m£ = rrtr or m£ = nie^fj,- 

In Table [5] we present numerical results for the normalized coefficient functions Ff{q‘^) 
and their averages. In columns 2 and 3 we list the values of F^{q^) for q^ = 50 GeV^ 
with zero and nonzero lepton masses. In order to avoid possible contamination from 
contributions of the '0 and T families we have chosen a q^ value in between these two 
families, namely q^ = 50 GeV^. This q^ value is small enough to highlight the helicity-fiip 
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i 

{nil = 0) 

{me = mr) 

W) i.'rni = 0) 

(Jf) {me = mr) 

1 

-0.9115 

-0.6257 

-0.3916 

-0.3491 

2 

-0.9115 

-0.9391 

-0.3916 

-0.3908 

3 

+0.9557 

+0.6561 

+0.6958 

+0.6537 

4 

+0.0030 

+0.0023 

+0.0203 

+0.0206 

5 

+0.0167 

+0.0132 

+0.0319 

+0.0319 

6 

+0.1875 

+0.1287 

+0.3589 

+0.3528 

7 

+0.0332 

+0.0228 

+0.2281 

+0.2284 


Table 2: Numerical results for the normalized coefficient functions at = 50 GeV^ 

and the average of P^{q^) over q^ G [dm^, [ttih — mzY] 

and lepton-mass effects in the vicinity of the threshold. On the other hand, this value of 
is far away enough from the threshold region where one would have to deal with the 
Coulomb singularity. We mention that the contribution of the 'ijj and T families to the q^ 
spectrum have been investigated in Ref. [10]. These contributions have been found to be 

small. 

Concerning the = 50 GeV^ values for the normalized coefficient functions, lepton- 
mass effects amount to —31 % for the functions ^ 13 ^ 7 , —21 % for the functions and 
-|-3 % for the function . The normalized coefficient functions quite small to 

start with. We mention that r-lepton mass effects are even larger for smaller values of g^. 

In columns 4 and 5 of Table [2] we also present average values {Tf) of the coefficient 
functions again for zero and nonzero lepton masses where the average is taken with regard 
to Y ■ In order to do the requisite g^ integrations one needs to include the relevant g^- 
dependent integration measure dehned by the differential (f distribution. Inserting the 
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necessary coupling and phase-space factors one obtains 


dV^ _ Bzii C^{q^) 9, 

dq^ d cos 9p d cos 6q dx ciw 8 ^ 


x-w^{q\ep,eq,x), 


where 


C^{q^) = 


9 


\pv(ml,q‘‘)\v. 


(17) 


(18) 


cos"^ 9w 4 • 15367r^ — m|)2 -|- m|r| 

Bzu is the branching ratio of the decay Z —)■ i.e. Bzu = r(Z —)■ i~^i~)/T{Z), where 

the rate for the decay Z —)■ {mip = 0) reads 

.2 1 


hz,, = r(Z ^ £+£-) = 


9 


-mz {vf + af). 


cos^ 9^J\l 1927r 

The magnitude of the momentum of the gauge bosons is given by 


\pv{p'‘,q^)\ = 


(19) 


( 20 ) 


where A(a, 6, c) = — 2ah — 2ac — 2bc is Kallen’s function. We then dehne partial 

differential rates according to 


dq"^ 


, BzaC^{ q^^ 

r'(p) 

^ew 


= 2 p^ 2q‘ 




( 21 ) 


The factors = m| and q^ are picked up when doing the integrations over the Z —)■ £'^£~ 
and Z* —)■ £^£~ phase spaces. The factor q^ is of crucial importance to cancel the 1/g^ 
singularity in the double spin-density matrix elements poo, Pot and pu at the lower end of 
the spectrum. The average values of the coefficient functions (Tf) can be calculated 
from the formula 




ldq'^2p‘^2q-^BzuC^(q-‘)Tpq^)IClP} ^ Ff 

r^' 


( 22 ) 


! dqnpnq^BzuC^{q‘)Ti{q^)ld£l 
The integration has to be done in the limits < g^ < {mn — ^zY- The denominator of 
Eq. fl2^ is nothing but the total rate T^ including lepton mass effects. When calculating 
the average values according to Eq. (12^ one can disregard all constant factors in the weight 
function C'^(g^). 


13 












Numerically, one finds that the averaged coefficient function {J^^) is by far the largest 
one. Lepton-mass effects are largest for (.Ffg) and amount to —10.8% and —6.1%. 

Using the average values {^i), the g^-integrated angular decay distribution can be 
written as 


dV'^ 


:(! +. 0 ). 


(23) 


2=1 


d cos 9p d cos 6 q dx Stt 
O ne can make contact with the work of Ref. [H] by taking the zero-mass limit miq —)■ 0 
of Eq. lET]) (for f = 0), neglecting the Z width and omitting the factor Bzu- In fact, using 
the notation mz = mzimii and (f = one obtains 


dT^ 


9 


mn 




dq^ cos'^ 6 w 4 ■ 30727r3 
in agreement with Refs. gmig. At = 0 one has 




(g2 - ml) 


dV^ 


if = 0 ) = 


9 


rriH 


C(i) (1 

^ew \ 


m 


2 \3 


= 4.02 ■ 10"^ MeV. 


(25) 


dg 2 cos^ dw 4 • 30727r^ 

Integrating Eq. fl24)) within the limits 0 < g^ < (1 — mzf one obtains the well-known 
expression nulls] 


r(i7 ^ z£+r) = 


9 


mn 


cos^ 6w 4 • 30727r^ 


F{mz) 


(26) 


with 


Fifhz) 


3(1 — Smi -|- 20mi) /3ml — 1 

-;- K -r- 77 ^—- arccos (-^— 

(4m| — 1)^/2 y 2m| 

- (1-m|)(^m| - ^-7 -3(1-6m| + 4m|)lnmz. (27) 

^ Zi Zi z' 


3.2 Single-angle decay distributions 

Integrating Eq. flTS]) over cos 9p and y and using Table HJ one obtains 

Wfql9q) = l(l + F^P2icos9q)) . (28) 

We define a convexity parameter Cf^\q‘^) as the second derivative of Eq. (l28jl with respect 
to cos 6q or, equivalently, as two times the coefficient of the cos^ 9q term in Eq. (|28|) . One 
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obtains 




3 fi + egf 3 Pu- 2pL 


{q^ - Am^,){v^ + aj) 


(29) 


2 ^ 2, fo + epQ 4 pu+L Ciwq^ + 2m}{yj+ ?,a}FgpslP u+l) 

In Fig. [2] we show a plot of the distribntion of the convexity parameter for both the 
mi = 0 and mi = m,- cases. Dne to the overall factor — Amj) one has C'f\q^) —)■ 0 
at threshold (maximal recoil) q^ = Amj, i.e. the cos^ 6 q distribntion is flat at threshold. 
This is clearly visible in the i = t case in Fig. O For i = e,p the vanishing of the 
convexity parameter C‘f\q^) at threshold is not discernible at the scale of Fig. |2l Instead 
Cj^\q‘^) —3/2 at threshold as mj —)■ 0 dne to the limiting behavionr of the two factors 

in Eq. fl2^ . The hrst factor goes to —3/2 becanse of the dominance of pi and the second 
factor goes to 1 in this limit. A closer look at the second factor in Eq. fl2^ reveals that it 
shows a steplike behavionr at threshold for mi —)■ 0 jnmping from 0 to 1. In fact, setting 
F| ps/pu+L ~ 1 in the second factor in Eq. (129|1 one obtains the limiting form 


(g^ - Amj){vj + a 




{q^ - Am])/{q^ + 2m|) . 


(30) 


Cewq"^ + 2 mj{vi + 3ajF^ps/pu+L) 

Eqnation flHHl) shows the advertised steplike behavionr as mj —)■ 0. Even for the mnon the 

deviation from 1 at = 10 GeV^ is a tiny one (0.67%). At the other end of the spectrnm 
at minimal (zero) recoil where g^ = {mu — the convexity parameter goes to zero in 
both cases since pf/ — 2pi ~ and = 0 at zero recoil. 

Fignre[2]shows that the convexity parameter is negative, i.e. the polar angle distri¬ 
bntion is described by a downward-open parabola which has a maximnm at cos Og = 0 with 
the maximal valne IF'^(g^, 6 q = 7r/2) = |(1 — C'f\q^)/3). In Fig. [3] we show a plot of the 
cosOq distribntion for g^ = 50GeV^. The cosOq distribntion is symmetric in cosOq dne to 
the absence of a term linear in cos 6q in Eq. (l28p . i.e. the distribntion is forward-backward 
symmetric. As expected from Eq. fl30l) the mi = m^ cnrve is considerably flatter than the 
= 0 curve. At cos 9q = ±1 the = 0 curve is close to zero since P2{cos 9q) = 1 at 
these points and J7f = fi/f^ ~ —1 due to the dominance of the longitudinal contribution 

pL- 
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Figure 2: distribution of the convexity parameter for me = 0 (dotted line) and me = m-r 

(solid line) for the process H ^ Z + Z*(—>- £+£“) 



Figure 3: cos0g dependence of the normalized decay distribution W^(q^,0q) for q^ = 
50 GeV^ in case of = 0 (dotted line) and me = m^ (solid line) 
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Figure 4: cos6*p dependence of the normalized decay distribution W^{q‘^,6p) for = 
50 GeV^ in case of = 0 (dotted line) and (solid line) 

Next we discuss the single-angle cos 6p distribution. From Table [U one reads off 

W^iq^Op) = ^ (l + ^iP2{cos9p)) . (31) 

The corresponding convexity factor is now given by 

_ 3 _ 3 fi + eg^ _ 3 pu - 2pL + 2m}{vj - QajFlps/P u- 2 l) , 22 ) 

^ 2 ^ 2 fQ + egQ 4 pu+L Ciwq"^ + 2mj{v‘^ + 3ajFgps/pu+L) 

The threshold value of the convexity parameter can now be seen to be given by C^^\q‘^) = 

—3/2 in both the = 0 and rrii = 771 ^- cases. We do not provide a plot of the convexity 

parameter C'f\q^) because lepton-mass effects are small even close to threshold. In Fig. |4] 

we plot the cos 6p dependence of the normalized single-angle distribution again for = 

50GeV^. There is practically no lepton-mass dependence in the cos 9p distribution. Since 

F 2 = Fi ill the zero mass case, the zero lepton-mass distributions in Figs. [3] and H] are 

identical to each other. 

Finally we turn to the normalized single-angle azimuthal distribution, where 










X 

Figure 5: x dependence of the normalized decay distribution x) for (f = 50 GeV^ 

in case of = 0 (dotted line) and = nir (solid line) 

In Fig. |5] we show the y dependence of W^{q^,x) again for = 50GeV^. The nonflip 
contribution to the coefficient of the cos 2% term clearly dominates the decay distribution 
since the leading contribution is given by ff ~ ajv‘^p+_ = ajv"^. The dominance of the 
cos 2x term is clearly evident in Fig. |5l Lepton-mass effects are generally small and amount 
to maximally ~ 1.2% at y = 0, 7r/2, 3/27r, 27r where the mass dependence mainly results 
from the normalization. 


3.3 The polarization of the off-shell gauge boson Z* 


In Sec. 13.21 we have already considered the single-angle cos 6 q distribution which we wrote 
in the form ~ (1 -|- P 2 {cos 6 q)). In this subsection we want to write the same angular 

decay distribution in terms of the transverse, longitudinal and scalar components pu, pL 
and ps of the double spin-density matrix Pmm'- One obtains 




= 2p^2q^BzuC^{q^) 


dq^d cos 9q 


3 3 

-(1 COS^ Oq) pu + x ^9 PL 

o 4 


(^(<?) 
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2 m^ 


- sin eqpu + - cos Oq Pl]v^ + -Fs{q )psa^ 


|(1 + cos^ eq)pu + j sin^ dq pl - 4ev^) 


= 2p^2q^BzuC^{q^) 

+3e(^Pu+LVe + Fl{q^)psafj 


( 34 ) 


Integrating the differential rate fl34)) with respect to cos 6 q, one obtains 

^ = 2pWB,aC^{P)\^ (ci«> + (35) 

Equation fl3^ can be seen to be the equivalent of the i = 0 piece of Eq. fl2l|) . Accordingly 
we define partial rates by writing 

= 2p^2q‘BzuC^(q'‘) (cW + pu.L , 

^ = 2p^2q^Bz,eC^g^)^Fl{g^)a-‘,ps. (36) 

In Eig. | 6 ]we display the dependence of the three partial rates dV^/dq^ (a = U, L, S) for 
the two cases rrii = 0 and mi = Lepton mass effects are largest for q^ values in the 
vicinity of the threshold. There is a substantial tauonic scalar rate for low values which 
partially compensates for the loss of longitudinal rate in the low g^ region. In the massless 
case, where there is no scalar partial rate, the longitudinal rate dominates the transverse 
rate up to ~ 700 GeV^. The transverse rate dVfj/dq^ shows a very small lepton-mass 
dependence. The transverse rates vanish at threshold due to their respective threshold 
factors. The g^ = 0 mass-zero rate can be seen to be in numerical agreement with the 
corresponding g^ = 0 value of Ref. [TT] listed in Eq. (I2H]) . 

The angular decay distributions are given by the normalized partial rates dVu^L^s/dq^ 
divided by the total rate dT[ 7 +^+ 5 /dg^. For brevity we denote these normalized rates by f7, 
L and S. Again we take our reference value g^ = 50 GeV^. The transverse-longitudinal- 
scalar helicity fractions are given by 


U : L : S = 0.06 : 0.94 : 0 U : L : S = 0.04 : 0.65 : 0.31. (37) 
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Figure 6: Differential rates dV^/dq^ (indices a = U, L, S and vanishing index for a = U+L) 
for the decay H Z{-^ e^e~) + Z*{-^ t^i~) with me = 0 and me = m^ 

The left and right three values refer to the modes H —)■ 2’(—)■ + Z*(—and 

H —)■ Z{-^ e+e“) + Z*(—r+r“), respectively. In the r mode one observes a substantial 
loss in the longitudinal rate which is compensated for by the appearance of the scalar rate. 
This has consequences for the cosOq distribution as shown in Fig. |3l 

In the upper part of Table [3] we list the total rate and the mean values of the transverse, 
longitudinal and scalar partial rates for the mass-zero modes and the r mode where the 
mean is taken with regard to q^. Lepton-mass effects reduce the total rate by 3.97%. 
The rate reduction is largest for the longitudinal rate where the rate reduction amounts to 
6.99 %. This is partially made up for by the appearance of the scalar rate which amounts 
to 4.14%. The average transverse rate is practically unaffected by lepton mass effects. 

3.4 Off-shell — off-shell decays H —> £+£“) + Z*(—> 

To conclude the section about Higgs boson decays into a pair of Z bosons, we briefly 
consider the case where both Z bosons are off-shell. The double-differential decay rate for 
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H ^ Z*{^ i+i-) Z*{^ i+i-) reads (4 ^ 4) 


dr^ 


dp^dq'^ 

X 


{p\q') = 


9^ 


8 • 192 ■ 1927r^ cos® 9iy rn^j 
1 1 


TTX^ . o o, I 

— \Vv{V )\VpVq 


(p2 _ 77^1)2 _|_ 777,|r| (g2 — m|)^ + m|r| 16 
X |Lf(p4Pr(p) + 3)F|(/)Lf(p4Po"^(p)| 


X Lf{q^)Pr{q) + 3P|4^)Lo4g4Po"4g) , 


(38) 


where we have written the result in terms of the spin-1 and spin-0 projections of the 


neutral current lepton tensor listed in Appendix B. The spin-1 and spin-0 propagators Pf^ 
and Pq'' are defined in Eq. (jl]). The velocity-type parameters Vp and Vq are defined by 
Vp = 2|p£p|/vF = and Vq = ‘l\viq\l = ^1 -4mfg/gT 

In writing down Eq. fl38|) we have chosen a p •H- g symmetric representation. This sym¬ 
metric form is very useful when one discusses the identical-particle decay H 
where two of the four contributing diagrams have the factorizing form of Eq. (138|1 . To 
achieve the p ■(-)■ g symmetry one has to add the scalar pieces to the p-side propagators in 
Eq. ([2]), i.e. one replaces Pf^(p) by Po®i(p) etc. The representation fl38|) is then obtained 
by expanding e.g. f dQqPZi(q) L\^Kq) Poiiiq) along Pf^ and Pq^. 

The contractions of the propagator factors can be calculated to be 


Prip)Pipu{q) := Pu+L{p\q') = 2 + 


2 ^ o , 0 , 1 , ^h\P2 


p2q2 


Prip)Popuiq) := psip^qy = -1 + 


{pqy 

p2q2 


= 2 + 1 + 

^hIpz^ 


p2q2 


p2q2 


Pr{p)P^M ■■= Pss{p\qy = 


{pqf 

p2q2 


= 1 + 


'm]j\pzA 

p2q2 


( 39 ) 

( 40 ) 

( 41 ) 


The transverse and longitudinal pieces in Eq. (I39|) are given by p;/ = 2 and pi = (pg)^/p^g^. 
The scalar-scalar contribution pss = Pq^{ p)Popu{q) in Eq. 1155]) appears multiplied by the 
product of helicity-flip factors mj/■ mj/q‘^ and can be neglected for all practical purposes. 
In the present case we take mip = rrie on the p side and meq = rur on the q side such 
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rz 

r§/r^ 

rf/r^ 

rf/r^ 

H ^ Z{^ 

e+e-) + Z*{^ i+i-) 




{mi 

= rrip) 

1.008 X 10-^ GeV 

0.4056 

0.5940 

0.0004 

{mi 

= mr) 

0.968 X 10-^ GeV 

0.4062 

0.5525 

0.0414 

H Z*{^ e+e-) + Z*{^ £+£- 

) 



{mi 

= rrip) 

2.449 X 10-^ GeV 

0.3879 

0.6119 

0.0002 

{mi 

= mr) 

2.405 X 10"^ GeV 

0.3881 

0.5936 

0.0183 


Table 3: Total and normalized partial decay rates for the four-body decays H —)• 
e+e“) -I- Z*{^ £+£“) (upper part) and H —)■ Z*{^ e+e“) -f Z*{^ i^i~) (lower part) 


that the symmetric appearance of Eq. (]38il is lost. In particular, we set Vp = 1 and take 
the mi ^ 0 limit of the hrst curly bracket in Eq. (I38ll replacing it by {vj + aj)Pi'^{p). 

In the lower part of Table [3] we present our numerical results for the off-shell - off-shell 
case. We list the total exclusive decay rate T^ and the averages of the partial decay rates 
Vfj/T^, rf/T^ and Tf/T^ for mi = m^ (hrst line) and mi = m-r (second line) on the q 
side {mi = miq). For the p side we specify to mip = mg. The off-shell - off-shell rates can 
be seen to be approximately twice as big as the on-shell - off-shell rates. The reason is 
that in the off-shell - off-shell case one picks up contributions from the peaking regions on 
both the p and q side. The helicity fractions remain practically unchanged except for the 
scalar contribution which is reduced by ~ 50 %. The reason is that the scalar contribution 
comes only from the q side whereas the normalizing rate is approximately doubled. Our 
result agrees with the result T^ = 1.0256 ■ 10“® GeV of Ref. |13] within 1.5 %. 

One can undo the smearing in Eq. fl38|) by the zero-width substitution 


TT 


(p2 — my)" 


+ mlrTl, 


-)■ 


myTy 


5{p^ - my) 


(42) 


22 








with V = Z. One then obtains 


dV^ 


9 


\Pv{P^,q^)\Vq 


X 


4 • ISShvr^ cos"^ 9w 
Bzu 


m 


H 


{Lf{q^)pu+L{ml, q^) + 3L^ {q^)Fl{q^)ps{m\, q^)] 


(g2 — rrt|)^ + m|r| 

BzamlC^(q‘^) {Lf(q‘)pu+L(’nz,q‘) + 3Lj'((/^)f|(j'^)ps(m|,9^)} , 


( 43 ) 


where dV^/dq^ denotes the differential rate into the g-side leptonic mode. As expected, 
this result coincides with Eq. fl35|) . The result for the (semi-inclusive) three-body decay 
H ^ Z + Z*(—)■ can be easily obtained from Eq. flT3l) by summing over all channels, 

i.e. by skipping the branching ratio factor Bz(.i- The result flT5]) without the factor Bzu 
can be seen to coincide with the result of Ref. [TT] . 

Finally, also the g-side secondary decay process can be considered to be exclusive. Using 
the decay rate for the decay Z —)■ in Eq. (IT^ for both the p and g sides, from Eq. (155]) 

one obtains 

d^pq , 2 23 ^ 1 g^m| \pv{p^,q^)\vpVq 

dp'^dq^ ^ 2 Sncos'^Ow + ofY 

BzuptTLz^ z BzuqmzT z p^q^ 

TT ((p2 — m|)2 -|- m|r|) TT ((g2 — m|)2 -|- m|r|) 

X |„J(1 + 2^)pr(p) + ajvlPr(p) + 3a] ■ 2^Fj(p2)P„'“'(p)| 

{ 777 ^ 777 ^ ') 

v](l + 2^)P,,,(q) + dlv]P,,M) + 3a] ■ 2-^F|(«2)P„„.,(«)|, (44) 

T^ denotes the rate into the exclusive leptonic modes on the p and g sides, respectively. 
In the product BzeipBzeeq there are terms diagonal and nondiagonal in flavour. The non¬ 
diagonal terms appear in pairs referring to the same exclusive channel. Thus one has to 
divide by a factor of two in Eq. fl4T|) as concerns the nondiagonal terms. The diagonal 
terms appear only once in the product. In the approximation that the interference contri¬ 
butions of the diagonal terms can be neglected, the factor 1/2 correctly counts the number 
of diagonal terms (see the discussion in Ref. ra). From the inclusive point of view the 
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Figure 7: Definition of the momenta p and q, the polar angles 9p and 9q, and the azimuthal 
angle x in the cascade decay H —)■ 

factor 1/2 in Eq. appropriately accounts for the identical particle factor of 1/2 in the 
rate. 

It is noteworthy that Eq. cannot be obtained in any gauge without considering 
the coupling of the off-shell gauge bosons to fermion pairs. If one calculates the rate for 
H —)■ Z*Z* in the spin-1 Lorenz gauge (also called Landau gauge), the result has to be 
multiplied by the effective factors p^/m| and q^lm\ in order to obtain Eq. (jllj). 

4 The four-body decay H —)► W~{-^ + 

Even though the yield of the {H iuiu) mode from Higgs decay is about 40 times larger 
than the yield of the {H —)■ mode, the identihcation of the H mode is 

much more difficult experimentally but can nevertheless be done [ 6 ]. In this section we 
write down the three-fold angular decay distribution of the cascade decay H IF“(—)■ 
^~V() -1- T^Vr)- In the charged-current decays there are no identical-particle effects 

such that one can also e.g. consider the decay H —)■ IF“(—)■ T~Vr) + As 

in the neutral-current case the r mass can be safely neglected on the on-shell side since 
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the scale is set by the W mass. At the end of this section we shall also discuss off-shell - 
off-shell decays where the r mass can no longer be neglected on the p side. A new feature 
appearing in the charged-current decays is the presence of a scalar-longitudinal interference 
effect which is parity-conserving but can mimic a parity-violating contribution. One can 
anticipate without explicit calculation that lepton-mass effects are not as important in 
the charged-current case since the corresponding helicity-flip contributions are four times 
weaker than in the neutral-current case. It is for this reason that we do not discuss the 
charged-current case in as much detail as the neutral-current case. 


4.1 Three-fold angular decay distribution for the 

four-body decay H —>■ W~{-^ + lT+*(—> t^Pt) 

Let us hrst consider the three-fold angular decay distribution of the decay H —)■ W~{-^ 

The polar angles 6 p and 6 q are dehned in the respective lepton pair 
center-of-mass systems while the azimutal angle x again describes the relative orientation 
of the planes as shown in Fig. [71 We use the zero lepton-mass approximation for the 
on-shell decay W~ —)■ but keep the r mass hnite for the off-shell decay W^* —)■ T^v-r. 

Again we split the angular decay distribution into its helicity-nonflip and helicity-flip 
part. Accordingly we write 

+ l'l-'v(9p.9,.x). (45) 

The nonflip decay distribution is given by 


{(P++ + P—)[l^ 

— 2(p++ — p _ )(cOS 6p — cos 9q){l — cos 9p cos 9q) 

— 4(p+o + P-o) sin 9p sin 9g{l — cos 9p cos 9g) cos y -f 2p+_ sin"' 9p 
+ 4(p_|_o — p-o) sin 9p sin 9q{cos 9p — cos 9q) cos y| 


9p sin^ 9q cos 2y 
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= m^). For the flip contribution one obtains 

9 

771 c 

+ 4poo sin^ dp cos^ dq — 2(p_|_+ — p _) cos dp sin^ dq 

— SFsPot sin^ dp cos dg + 4:FgPu sin^ dp — (p+o + p_o) sin 2dp sin 2dq cos x 
+ 2(p+o - P-o) sin dp sin 2dq cos x + ‘^Fs{p+t + P-t) sin 2dp sin dq cos x 
— AFs{p+t — P-t) sin dp sin dq cos x — 2p+_ sin^ dp sin^ dq cos 2x|, (47) 

including the extra minus sign for the spin-1 - spin-0 interference contributions linear in 
Fs- We have used the velocity-type parameter Vq = 2,\'piq\l = 1 — mjq/q"^. 

In Eqs. fl46l) and fl47|) we have also included the contributions from the parity-violating 

terms proportional to (p++ — p _ ), (p+o — P-o) and (p+t — p_t)- These coefficient functions 

are not populated by the parity-conserving SM (HVV) coupling. In Appendix A we briefly 
discuss the contribution of a parity-violating non-SM coupling proportional to 
which would populate the (p++—p_), (p+o —P-o) and {p+t — P-t) coefficient functions [T71 


Again we write the result in terms of the Legendre polynomials Pi (cos 0) = cosd and 
P 2 (cos 6^) = I (3 cos^ — 1), 


{2p^2qY^W^{0p,dq,x) 


IQVq 

9 


10 


i=0 


ISUg 

9 


J:uT' + e9r)h.(e,,e„x). 


i=0 


(48) 

The coefficient functions and can be found in Table |4] where we have now dropped 

the non-SM contributions proportional to (p++ — p_), (p+o — P-o) and (p+t — p-t)- It is 

noteworthy that three new angular structures proportional to cos dq {i = 8, 9) and sin dq 
{i = 10) are generated by a helicity-flip contribution. The hrst of these contributions 
{i = 8) give rise to a nonvanishing forward-backward asymmetry in the cos dq distribution 
as discussed later on. 
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i 

fY 

W 

gY 

hi{9pi 9q^ x) 

0 

PU+L 

\pU+L + iFjps 

1 

1 

\{Pu - ‘^Pl) 

-\{pu - ‘^Pl) 

P 2 (cOS 9q) 

2 

\{pu - ‘^Pl) 

\{pu - 2 pl) -\FIps 

P2{cos 9p) 

3 

\{Pu + ^Pl) 

-\{Pu + ^Pl) 

P2{cOs9p)P2{cOs9q) 

4 

-4 Pc/ 

0 

cos 9p cos 9q 

5 

“|(p+o + P-o) 

0 

sin 9p sin 9q cos y 

6 

^(p +0 + P-o) 

“^(P +0 + P-o) 

sin 29p sin 29q cos y 

7 

ip+- 

-|P+- 

sin^ 9p sin^ 9q cos 2 y 

8 

0 

—3Fs Pot 

cos 9q 

9 

0 

3-^5 Pot 

P 2 {cOS 9p) cos 9q 

10 

0 

1 Fs{p+t + P-t) 

sin 29p sin 9q cos y 


Table 4: Coefficient functions appearing in the three-fold angular decay distribution of the 
decay H i-p^) + W+*{^ t+v^)- 


We dehne a normalized decay distribution 


w^{e,,e,,x) = 


jw^{e’p,e’^,x')dcose’pdcose’^dx' stt 

where (and fT = ~9Y = and where 

= fY + ^qY = Pu+L + \e{pu+L + 3F|p5). 

The differential decay rate distribution is given by 

dV^ „ C^(g2) 


(l+i:-Fffe,(»„»„x)). (49) 


(50) 


dq^ d cos 9p d cos 9g dx 


= B 


Wiu- 




Stt 


16n, 


(51) 
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where 


C^{q^) = 


g\ \pv{mw,q^)\vq 


(52) 


15367r3 rriH {q^ - 

(note the additional factor Vq in the numerator) and The decay rate 

for W ^ {rrii = 0) is given by 

^2 


'^wtu = r(iy ^ ^v) = 

4o7r 


(53) 


Partial differential rates are dehned according to 


^ 2q^Bwe.C^{q^W{q")- (54) 

The average values of the coefficient functions is given by 

l<lqV2q‘B„,,C"‘\q‘)Tmq‘) rr 
' ■ ' !dqnpnq^Bw,,C'^{q^)F”{q^) T'*'' 


where the integration over q^ runs from to {mn — mwY- 

In Table |5] we present numerical results for the normalized coefficient functions (g^) 
and their averages. In columns 2 and 3 we list the values of (g^) for g^ = 50 GeV^ 
with zero and nonzero lepton masses. Lepton-mass effects amount to —16% for the func¬ 
tions g —11.8% for the functions and only -|-0.4% for the function . The 
normalized coefficient functions are zero for zero lepton masses. For mi = m^- the 

coefficient functions become quite large at 16.1%. Again, lepton-mass effects would be 
even larger for smaller values of g^. Compared to the H ^ ZZ* case, lepton mass effects 
are smaller by approximately a factor of 1/2. 

In columns 4 and 5 of Table [5] we also present average values ) of the coefficient 
functions again for zero and nonzero lepton masses where the average is taken with regard 
to g^. On average lepton mass effects can be seen to be quite small. 






i 

= 0) 


(T-f) (m, = 0) 

(Jf") {me = mr) 

1 

-0.9549 

-0.7983 

-0.3965 

-0.3817 

2 

-0.9549 

-0.9585 

-0.3965 

-0.3985 

3 

+0.9775 

+0.8172 

+0.6983 

+0.6809 

4 

-0.6759 

-0.0603 

-0.9052 

-0.9006 

5 

-0.5431 

-0.4847 

-1.4306 

-1.4205 

6 

+0.1358 

+0.1135 

+0.3576 

+0.3533 

7 

+0.0169 

+0.0141 

+0.2263 

+0.2244 

8 

-0.0000 

-0.1614 

-0.0000 

-0.0176 

9 

+0.0000 

+0.1614 

+0.0000 

+0.0176 

10 

+0.0000 

+0.0015 

+0.0000 

+0.0029 


Table 5: Numerical results for the normalized coefficient functions (g^) at = 50 GeV^ 
and the average of (g^) over g^ G [mj, {mn — ttiwY] 
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4.2 Forward-backward asymmetry of the lepton 

Let us take a closer look at the cos 6q distribution determined by the coefficient functions 
and ■ The normalized cosQq distribution is given by (see Table 01) 


Oq) = ^ (l + J^rP2{cOSeq) + COsOq) 


(56) 


Contrary to the H ^ ZZ* case one now has a contribution linear in cos 9q which implies 
a nonvanishing forward-backward asymmetry. It is interesting to note that the source of 
this parity-odd term in Eq. flFr|) results from a parity-conserving interaction. Consider 
the content of the currents coupling to the W*-. and 0“). The 

scalar-longitudinal interference contribution leading to a nonvanishing forward-backward 
asymmetry can be seen to result from the parity-conserving interference of the products of 
currents l/(0+)F(l“), 4(0“)4(1+). We mention that a parity-violating {HW~^W~) cou¬ 
pling as discussed in Appendix A would also give rise to a nonvanishing forward-backward 

asymmetry proportional to (p++ — p _) (see Eq. f06l) L 

In Fig. Owe display the normalized cos 6q distribution for a hxed value of = 50 GeV^. 
Since pu ^ Pl for = 50 GeV^, the governing feature of the distribution is described by 
a downward open parabola. The convexity parameter is proportional to (1 — s), leading 
to a smaller convexity for the mi = rrij. distribution as can be seen in Fig. 01 There is 
a pronounced forward-backward asymmetry in the r mode. According to Eq. fl56l) the 
forward-backward asymmetry of the cosOq distribution is given by 


4 ^f-Tb mj ^ -3Fs{q^)pot{q^ 


(57) 


F L B ^ q 

Since pot and Fs are positive, the forward-backward asymmetry is negative as also 
shows up in Fig. 01 In fact, one calculates 


ApBiq"^ = 50GeV2) = -0.081. 


(58) 


For smaller values the forward-backward asymmetry becomes even more pronounced. 
In the last column of Table El we list the average value oi Afb which is given by (Afb) = 
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Figure 8: cos 9 g dependence of the normalized decay distribution {q‘^, 9 g) for = 
50 GeV^ for = 0 (dotted line) and mi = 171^ (solid line) 


average, the forward-backward asymmetry is much smaller than for = 

SOGeV^. 


4.3 The polarization of the off-shell 


Let us rewrite the single-angle cos 9q distribution in Eq. fl56l) in terms of the transverse, 
longitudinal, scalar and the scalar-longitudinal interference contributions. One has 


dr 


w 


dq^d cos 9q 


+ 


2g2 


= 2p^2q^Bwe.uC^{q^)W{l + cos^ 9q)pu + ^ sin^ 9qPL 


- sin^ 9qpu + - cos^ 9qPL - ^Fs{q^) cos9qpto + -Fg{q'^)ps 


Integrating the differential rate (IMD with respect to cos 9q, one obtains 

= 2p^ 2q^ BwiuC^iq^) |(P(7 + Pl) + ^ {pu + Pl) + ^Fg{q^)ps 
Partial decay rates are accordingly defined by 


dq'^ 


= 2p^ 2q^ BwiuC^iq^) ^ j Pu,l , 


(59) 


(60) 
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0 500 1000 1500 2000 

q' [GeV'l 

Figure 9: Differential rates dV'^ / dq^ (indices a = U, L, S and vanishing index for a = U+L) 
for the decay H W~{-^ e~Ue) + with mi = 0 and mi = m^ 

^ = 2p^2fB„.C«'(q^0Fl{f)ps. ( 61 ) 

In Fig. [9] we display the q"^ dependence of the partial decay rates. Lepton-mass effects 
show up mainly for lower values between threshold and approximately 150 GeV^. For 
example, at = 50 GeV^ the helicity fractions of the 1F+* change according to 

U : L : S = 0.030 : 0.970 : 0 ^ U : L : S = 0.028 : 0.893 : 0.079 (62) 

when going from the e,/i modes to the r mode. The picture is similar to the H ^ ZZ* 
case. However, lepton-mass effects are less pronounced in the H —)■ WW* case. 

In the upper part of Table Owe list the total rate and the mean values of the transverse, 
longitudinal and scalar partial rates for the mass-zero modes and the r mode where the 
mean is taken with regard to q^. In this case, lepton-mass effects reduce the total rate 
by 0.76 %. The rate reduction is largest for the longitudinal rate where the rate reduction 
amounts to 2.1%. On average, the rate reduction for (F^) = F^^/F is still a considerable 
1.3 % while the average of the transverse rate is practically unchanged. As in the W ^ ZZ* 
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case, the loss of longitudinal rate is mainly compensated for by the appearance of the scalar 
rate. Our result agrees with the result = 2.4135 • 10“’^GeV of Ref. [13] within 1.6%. 


4.4 Off-shell — off-shell decays H ^ W i ui) + 


Again, we conclude the section by considering the case where both W bosons are off-shell. 
Using again the narrow-width approximation (1421) for V = W also on the p side, the 
exclusive off-shell - off-shell rate is obtained from the double integral 


-pW 

^ pq 




dp^ 


Bwiup w 


rirriH-y/^)^ 


7r((p2 - 


£q 


dq^ 


Bwiuq w 


7r((g2 - 


■^w 
0 > 

(63) 


where denotes the rate into the exclusive leptonic modes on the p and q sides. The 
total inclusive mode is obtained by summing over all exclusive modes including the quark- 
antiquark modes, i.e. by setting Bwivp = Bwtuq = 1- The rate function T^ in Eq. (I64]l 
reads 


0 


VpVq 


^\pv{p^,q^)\— -2 a A 

Svr m^mjj 64 


Lr(p2)pr(p) +3Pi(p^)C(p^)Po"^(p) 


X Lr{q^)Pr{q) + ^Fl{q^)L'^{q^)Pr{q) , 


(64) 


where Vp = 1 — rrilp/p^ and = 1 — mjg/q'^. and are the spin-1 and spin-0 
projections of the charged current lepton tensors listed in Appendix C. When the lepton 
masses are taken to be zero, the rate function (l64)) simplihes to 

{me = 0) = ^\pw\-4^ Pu+L (65) 

oTT TTt-^TTl ^ 

which agrees with the results in Refs. naE]. 

In Table [6] we have listed our numerical results for the off-shell - off-shell rates and the 
averages of the polarization of the gauge bosons. The off-shell - off-shell rates approxi¬ 
mately amount to two times the off-shell - on-shell rates. The reason is again that one 
picks up contributions from the peaking regions both on the p side and on the q side. Put 


33 










^w/^w 

pW j^W 

pW/pW 

(^fb) 

H W-{- 

-^e-pe) + W+*(-^ i+ue) 




{me 

= rrif,) 

4.754 X 10-®GeV 

0.4023 

0.5976 

0.0001 

-5 X 10-5 

{me 

= m^) 

4.718 X 10-®GeV 

0.4033 

0.5894 

0.0073 

-0.0059 

H W-*{ 

e-pe) + W+*(^ 

i+ue) 




{me 

= rrif,) 

1.009 X 10-5 GeV 

0.3952 

0.6048 

3 X 10-5 

-2 X 10-5 

{rrii 

= rrir) 

1.005 X 10-5 GeV 

0.3956 

0.6009 

0.0035 

-0.0029 

H W-*{ 

i-Pe) + W+*(^ i+ue) 




(me 

= "i/.) 

1.009 X 10-5 GeV 

0.3952 

0.6048 

6 X 10-5 

-6 X 10-5 

(me 

= rrir) 

1.001 X 10-5 GeV 

0.3960 

0.5970 

0.0070 

-0.0076 


Table 6: Total and normalized partial decay rates and the average value of the forward- 
backward symmetry for the four-body decays H —)■ e~Pe) + (hrst 

part), H —)■ W~*{^ e“Pe)+hh+*(—)■ (second part), and H —)■ W~*{^ i~Pi)+W'^*{^ 
(third part) 

in a different language the off-shell - off-shell rate (|63ll corresponds to the sum of the two 
off-shell - on-shell rates H —)■ W~W~^* and H —)■ W~*W^. The polarizations listed in 
Table |6] are thus an average of the respective polarizations on the p side and the q side. 
This also explains the fact that the scalar polarization is only one-half of the on-shell - 
off-shell value since it is contributed to only by the q side H —)■ W~W^* channel. This is 
no longer true for the H —)■ hh“*(—>■ T~Pr) + t^^t) rnode where the scalar rate 

obtains contributions from both the p side and the q side. 

The off-shell - off-shell rates slightly exceed twice the off-shell - on-shell rates which 
provides a measure of the quality of taking the zero-width approximation on the on-shell 
sides. Compared to the rates for H Z*Z* listed in Table [3l lepton-mass effects can be 
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seen to be more than five times smaller in the H ^ W *W~^* case. 

5 Summary and conclusions 

We have discussed lepton-mass effects in the rate and the angular decay distributions in the 
four-body decays H —)■ Z{^ + Z*{-^ r+r“) and H —)■ W~{^ i~p£) + W^*{^ 

where the gauge bosons Z and W~ are on their mass shell. Lepton-mass effects are larger 
for the H —)■ ZZ* mode where we find a reduction of 3.97% in the r rate relative to the 
e, fi rates. In the H —)■ WW* case the rate reduction of the total rate is smaller. In this 
mode we find a rate reduction of 0.76% relative to the zero mass case. Differentially, the 
rate reduction through lepton-mass effects is significantly larger at the lower end of the 
spectrum in both cases. For both modes we find a significant reduction of the longitudinal 
rate through lepton-mass effects in the lower region from threshold to 200 GeV^. In 
this region the transverse-longitudinal composition of the off-shell gauge bosons is consid¬ 
erably changed. The reduction of the longitudinal rate in this region is partly compensated 
for by a significant scalar contribution. In the charged-current case one finds a nonvanish¬ 
ing forward-backward asymmetry in the cos 6q distribution through lepton-mass induced 
scalar-longitudinal interference effects. The forward-backward asymmetry can become 
quite large in the low-g^ region. 

We have also discussed the case when both gauge-bosons go off-shell. Double smearing 
with the appropriate Breit-Wigner functions increases the overall rate. Lepton-mass effects 
become weaker in the double smearing process. 

We have employed helicity methods in our analysis which has allowed us to present our 
analytical results for angular decay distributions and partial rates in compact form. In 
particular, the inclusion of lepton-mass effects in the helicity formalism is straightforward. 

Experimentally it will not be so simple to identify the r modes in the four-lepton 
decays of the Higgs. In this context we mention that the detection efficiency for r leptons 
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in their hadronic decay channels at the LHC is being continuously improved (see Refs. |l5l 
Eg ST]). Nevertheless, an accurate Monte Carlo event generator for decays involving the 
r leptons should include lepton-mass effects for which we have supplied the appropriate 
matrix elements in this paper. This would e.g. be relevant for modelling Z —)■ rr processes 
as background for the search for the decay H —)■ [HI ES [50] . 

Acknowledgments 

This work was supported by the Estonian Research Council under Grant No. IUT2-27. 
J.G.K. would like to thank A. Denner, B. Jager, M. Rauch, H. Spiesberger, B. Stech and 
J. Wang for useful discussions. S.G. acknowledges the support by the Mainz Institute of 
Theoretical Physics (MITP). 


A Helicity amplitudes for H —)► VV* 


The helicity amplitudes for the transition H —)■ VV* are dehned by 


U _ ^ 

^mn — ^ 


{Xv,p)Haa'£*“ (Av*,g), 


(Al) 


where s*°‘{Xv,p) and 8*°"'{Xy*, q) are the respective polarization four-vectors on the p side 
(on-shell side) and on the q side (off-shell side). We shall evaluate the helicity amplitudes 
in the Higgs rest frame with the 2 ; direction dehned by the direction of the off-shell V* 
boson. One therefore has to rotate the polarization four-vectors on the p side by 180° 
which we indicate by the “bar” symbol. 

The respective polarization four-vectors in the Higgs rest frame are given by 


on-shell side 
off-shell side 


^“(±,P) = ^(0;±l,-c0) £"(0,p) = — (|py|;0,0,-po) 

V 2 TTlv 

£^{±,q) = ^(0;Tl,-b0) ef^{0,q) = ^(|pv|; 0, 0, go) 

= ^(?o;0,o, Ipv-I) 
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where po = ~ <?o = {fn'ii + ~ p‘^)/i‘^^H) and 

^{pqY -p‘^q^ 


1 

2mH 


\Pv\ = 7 ; -= 


niH 


(A2) 


such that = (po; 0, 0, —\pv\), q^ = {qo', 0, 0, \pv\)- It is convenient to avail of the covariant 
representations of the longitudinal and scalar polarization four-vectors. They read 

1 / „ s T)^ 


on-shell side : £^(0) = 


{{pq)p^ - P^q^) ^ 




Vt^\J{ pqY - p^q^ 

off-shell side : £^(0) = ^ rr—.^ {{pq)q^ ~ q^P^^ ^ 


V¥\fipq)^^-^ 

The helicity amplitudes can then be calculated to be 




(A3) 


//++ = H__ = 1, Hoo = Hu = 


pq 




Hot = H,o = 


V jpqy - p'^q^ _ m^ipvi 

(A4) 


The coefficient functions are written in terms of bilinear forms of the helicity 

amplitudes for which we choose the following abbreviations 


Poo — l-Hool ) 


hfnnl^, p±± — |hf±±P; Ptt — ReTToi-^oo 


P±o 

= ReH±±HQQ, 

Pin 

= Reii±±ii;^, 


Pit 

= ReH±±HQu 

Pti = 

= RefiotR±±, 


Pot 

= Re if 00 -Rot) 

Pto = 

Re HioHqq. 

(A5) 


We sometimes refer to these bilinear forms as the double spin-density matrix elements of the 
gauge boson pair since the bilinear forms describe the entangled polarizations components 
of the gange boson pair. The SM values of the double density matrix pmm' are given by 


-1 . _ 

/^ib=F ? Poo 9 2 

/V-W^ 


P++ = P— = 1, 

pq ^{pqY - p'^q^ _ mH\pv\ 


^ ni] 


Pot — Pto — 

Ptt = 


p2q2 

^2 ^ 2^2 


2myq^ 

(pqf-p‘q‘ m% ^ 


. . = (i + -^IAf). 

\ q^rriy ^ 

(2 2 2\ 

[mu -ruy-q j , 


p2g2 


2 2-\Pv\^ 

q^rriy 


(A6) 
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One notes the following SM relations 


Ptt — Poo ~ 1) Pot — \fPmfht) P±o — P±t — \/^- (A7) 

At maximal recoil where and/or tend to zero, the dominant double spin-density 
matrix elements are poo = Pot = Pto = Ptt- At minimal recoil = {niH — rnyY where 

|pV| —)■ 0, the dominant contributions are p++ = p _= poo while poi and pu tend to zero. 

Some authors prefer to use Cartesian components for the transition matrix elements [H] 
instead of the helicity components used by us. The relation between the two representations 
is given by 

All = + Ao = Hoo- (AS) 

We see no particular advantages to write the angular coefficient functions in terms of their 
Cartesian components. 

Battacherjee et al. considered two additional non-SM (HVV) coupling structures [TS] . 
They write down the effective coupling structure 

V^^•' = - {pq)g^'') + ice^^^'P^pq,. (A9) 

The helicity components are then given by 

^00 = ~ H±± = a =F c^J{pqf - p'^q^. (AlO) 

There are no contributions of the new coupling structures to H^f It is clear that one now 

has a contribution to the difference (Tr++ — H _) resulting from the parity-violating term 

proportional to e^'^P'^Ppqa. 

B Helicity representation of the 
neutral-current lepton tensor 

We calculate the helicity representation of the neutral-current lepton tensor on the off-shell 
q side. The corresponding expressions for the on-shell p side can be obtained by setting 
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the lepton mass to zero and replacing q ^ p. We work in the center-of-mass system of 
the lepton pair with the z direction dehned by which we refer to as the helicity system. 
The kinematics in the helicity system is given by 





£'^(0) = (0;0 

where 

= 1 — -■ 

tVV 

^ fjLU 

= tv(7„(7- 

T AA 

= Tr(7p75(/’' 

tVA 

= i7 = Tr( 


= i. //rSC 
2 






The total neutral current lepton tensor is composed according to 

- 2veaeL]^^ + 

where the neutral current is dehned by with 

= —1 + 4 sin^ ai = —1 for i = e, p, t. 


(B2) 


(B3) 


(B4) 


In order to calculate the helicity representation of the lepton tensors in the helicity 
system one needs to evaluate 




(B5) 


All objects referring to the helicity system are denoted by a hat symbol. The contractions 
are done in the helicity system using the representations (IBlll . Using the explicit forms (IBip 
one calculates 


jVV{q) _ 

jVV(q) _ 
-^00 — 

-J AA{q) _ 


2g^ 


jAA{q) 


= 2g^(l-4e), 


jVA{q) _ fAV(q) 


= T2g'(l-4£)i/^ 


8q^e, 


YAA(q) _ PI 

-^00 ~ U, 


£7'” = 0, 


2^ fVV,AA{q) _ jVV,AA(q) _ p, 
n+n — — U- 


(B6) 
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The ratio of helicity-flip and helicity-nonflip contributions can be seen to be given by 
L57/LS = 4e and Lr/r±t = 4e/(l - 4a). 

The components of the off-shell lepton tensor Lx^x'^ in the Higgs decay system [z axis 
along the Z* direction) needed in the evaluation of Eq. ([8]) can be obtained by rotating the 
components of the lepton tensor y dehned in the helicity system according to (double 
indices are summed) 


n,Ai(9".cos«„x)= E (B7) 


<i°m = 1, <„.(«) = 


J=0,1 

where Wigner’s d functions are given by 

^ i(l-|-cos6*) —^^sin^ |(1 —cos6*) ^ 

^sin^ cos 6* —^sin6* 

^^|(l-cos6*) :^sin6' i(l-Fcos6*) y 

The rows and columns in the spin-1 part of fIBSp are labelled in the order (-1-1, 0, —1). 
One obtains 

= 4eal L'S = ig = 0, L\f = ig = 0, 


(B8) 


^2\-l T^q) _ 


(2g^)-L±± 


00 


= |(1 + cos^ d){vj + ajvq) ± 2v£aev cos 6 + 2ev‘j sin^ 9, 
= {vj -1- alvg) sin^ 9 + cos^ 9, 


{2q^) ± (n^-I-a^n^) sin26'— 4n£a£nsin6* =F 4£n|sin26*) 


{2q 


T(q) _ r Wt 




^±T 


(i(n| -h a^nj) sin^ 9 - 2evj sin^ 9) (B9) 

The corresponding expressions for the on-shell side lepton tensor Lp 2 Xz\'^{cos9p) can again 
be obtained by rotation. Note that in this case one has y = 0 as is evident from Fig. [1] 
The spin-1 and spin-0 projections of the neutral lepton tensor needed in the main text 
are given by 


Li(g2) = = Tf/+i(g2) = 4g2(^n|(l + 2^) + a^(l - 4^)), 

Lo(g^) = P^''{q^)L^^{q^) = Ltt{q^) = Aq'^ape, 


(BIO) 
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where the spin-1 and spin-0 projectors read 


^r(g^) = -9^^ + pr{9^) = (Bii) 

Similar relations hold for the p side. 

C Helicity representation of the 
charged-current lepton tensor 

The lepton tensors are given by 

on-shell side : Ljfj = Tr - 75)^7^(1 - 75 )) 

= 8 + \p^giiy - , 

off-shell side : Ljfj = Tr (/+ + - 75 )l^ 7 v(l - 75 )) 

= 8 + it\{<f - ^l)9pu + iepupaq'^i^'') ■ (Cl) 

The kinematics for the off-shell q side is given by 

g“ = y/^(1; 0, 0, 0) i+^ = + e-,0,0,1 - e) z/“ = - £)(1; 0, 0,-1) 

£''(0) = (0; 0, 0,1) £^(±) = ^(0; ^1, -b 0) e^{t) = (1; 0, 0, 0) (C2) 

(for the on-shell p side set e = 0 and q ^ p)- The nonvanishing components of the helicity 
representations of the lepton tensors can then be evaluated to be 

on-shell side : Ttlt = 8m^ 

off-shell side : L++ = Sq'^Vr = L^ot = — ^if = 4:'mlvr (C3) 

where Vr = 1 — rri^/q^ = 1 — e. Note that the ratio of helicity-flip and helicity-nonflip 
contributions are now given by e.g. LqoV^++ = rn^/2q^ = e/2. 

As in the neutral-current case the components of the lepton tensor in the Higgs 

decay system {z axis along the hh+* direction) needed in the evaluation of Eq. ([8]) can be 
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obtained by rotating the components of the lepton tensor in the helicity system 

according to (double indices are summed) 






(C4) 

'"VV ''VV ' 'lA/ ''\A/ '^'1/1/ >''Tir 

J,J'=0,1 

One obtains (the rows and columns of the matrix are ordered in the sequence {t, +1, 0, —1)) 


{2q\) ^Lx^y^{cos9,x) = 

/ n 0 


0 


sm 


0 

2 


0 

0 (1=fcos6')^ =F^(1 =F cos 6*) sin6*e*^ 

0 =F;^(1 =F cos6*) sin^e"*^ 2sin^6' =f^( 1 ± cos6*) sin6*e*^ 

=F^(1 ± cos6*) sin6*e“*^ (licos^)^ 




0 


sin^ 6e 


/ 


/ 


+e 


— ^sinOe'^^ 


2 cos 9 ^ sin 6*e*^ ^ 


— ^sm9e 
2 cos 9 


sin"^ 9 —h sm 


-^sin20e-*^ 


sin 26*e*^ — sin^ 

sin 26*e*^ 


2cos2 0 


(C5) 


/ 


y ;^sin6*e — sin^6*e ^sin26*e sin^6* 

The upper/lower signs refer to the decays W~ —)■ i~ui and W~^ —)■ The corresponding 

expressions for the on-shell-side lepton tensor y (cos^p) can again be obtained by 
rotation. However, in this case one has x = 0 as is evident from Fig. [71 Further, one has 
to use the lower signs in the matrices flC5F 

The spin-1 and spin-0 projections of the charged lepton tensor are given by 


LYU/) = = 8,^ (1 - T) (l + p) , 

lY('/) = a'“'(9tFD = 8,yi-T)ie, 


(C6) 
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